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The de Haas-van Alphen effect (dHvAe), describing oscillations of the magnetization as a function of mag¬ 
netic field, is commonly assumed to be a definite sign for the presence of a Fermi surface (FS). Indeed, the 
effect forms the basis of a well-established experimental procedure for accurately measuring FS topology and 
geometry of metallic systems, with parameters commonly extracted by fitting to the Lifshitz-Kosevich (LK) 
theory based on Fermi liquid theory. Here we show that, in contrast to this canonical situation, there can be 
quantum oscillations even for band insulators of certain types. We provide simple analytic formulas describing 
the temperature dependence of the quantum oscillations in this setting, showing strong deviations from LK the¬ 
ory. We draw connections to recent experiments and discuss how our results can be used in future experiments 
to accurately determine e.g. hybridization gaps in heavy fermion systems. 

PACS numbers: 


Introduction. Landau quantization of electrons [1], which 
leads to quantum oscillations (QO) of physical observables as 
a function of applied magnetic field [ 2 ], has been one of the 
cornerstones of condensed matter physics. On the one hand, 
it leads to new phenomena such as the integer quantum Hall 
effect [4] and its fractional version [5]. For the latter, it even 
induces an unexpected new phase of matter beyond the stan¬ 
dard Landau classification [ 6 ], which ignited the field of topo¬ 
logical phases [7]. On the other hand, it is itself an invaluable 
tool for the characterization of correlated metallic systems [3]. 
The canonical LK [ 8 ] theory of QO in metals showed that the 
periodicity, e.g. of the magnetization, is proportional to ex¬ 
tremal cross sectional areas of the FS, thus turning QO into 
a precise quantitative and by now standard tool for determin¬ 
ing FSs. In addition, Lifshitz and Kosevich showed that it is 
possible to study correlation effects by extracting the effec¬ 
tive mass, m*, from the temperature dependence of the QO 
amplitudes given by (for the first harmonic) 
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and the cyclotron frequency ujc = ■ 

Later the LK theory was extended to include more general 
self energy interaction effects [9-12], but these always pre¬ 
served the general structure of the LK theory only renormal¬ 
izing parameters, e.g. m*. It still comes as a great surprise that 
experimentally almost all materials, from weakly interacting 
metals to strongly correlated heavy fermion systems [13-15] 
or copper oxide high temperature superconductors [16-20], 
are consistent with a LK description which is manifestly an 
effective single particle theory. There have been only very 
few exceptions for heavy fermion systems, e.g. CeCoIns [21] 
and most recently the tentative topological Kondo insula¬ 
tor SmBg [ 22 ], violating the general temperature behaviour, 
Eq. (1). There have been recent theoretical studies on QO 
which explored novel effects due to symmetry breaking from 
commensurate [23, 24] or incommensurate [25] charge den¬ 
sity waves but they remained in the canonical LK framework. 
A notable exception is given by Ref.26 which derived a gener¬ 


alized formula for exotic quantum critical systems described 
via non-perturbative field theories. 

Historically, the firmly established understanding of QO is 
tied to the existence of a FS, which in principle impedes the 
following simple question: Can there be QO in an insulator? 
In this Letter we show that, surprisingly, the general answer 
is yes. This arises if the cyclotron frequency huJc is of the 
order of the electronic gap and the band structure picks out a 
particular area of the Brillouin zone (BZ), as described below. 
We further show that, even in this non-interacting setting, the 
electrons exhibit anomalous non-LK QOs. 

We show that a simple band insulator of itinerant electrons 
hybridized with a localized flat band does exhibit well-defined 
QO. The periodicity is given by the area defined by the inter¬ 
section of the unhybridized bands even if the chemical poten¬ 
tial, /i, is inside the hybridization gap or inside the flat part 
of the FS. In the latter case, the periodicity is equally unusual 
because it is not proportional to the FS area. We And that the 
temperature dependence of the oscillation amplitudes strongly 
differs from the standard LK theory: First, if p, is inside the 
gap QO amplitudes have a maximum at a temperature set by 
the hybridization gap, 7 . Second, for a chemical potential 
inside the bands but close to the flat regions the behaviour 
is even more complex and governed by an additional energy 
scale, 6p, which is the distance of p above the bottom of the 
upper band. For Sp < 2^ there is a characteristic steep 
increase of the amplitudes towards lowest temperatures. 

Our main result is the general temperature dependence 
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A simple approximate formula 
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is valid in the regime huic ^ 27 or more generally for T ^ 
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= 2 k “ f fu^c{l + 1 ) with Wc = eB/m, 

and ^ j: ^ 7V$ with 7V$ = the number of flux quanta 
$0 = 11 through the system area A. We have neglected the 
Zeeman energy splitting of spin components. For each LL 
index I we have two energies with E-{1) < E^{1) for all 1. 
Note that for the lower band E-{1 00 ) W, giving a 

divergent density of states; this is an artefact of the continuum 
flat band which needs to be regularized. 

Anomalous de Haas-van Alphen effect. We calculate the 

magnetization M from the grand canonical potential (A:^ = 1 ) 
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Figure 1: (color online.) Main figure (i), Quantum oscillations of the 
magnetization, M, as a function of VF/feuc oc 1/B. Inset (ii), sketch 
of the band structure for our model (exaggerated hybridization gap 
for better visibility) and positions of the different chemical potentials 
p. If p, is far away from the gap (black dashed and dot dashed), 
which is opened by hybridizing a localized flat band with an itinerant 
band, the periodicity of standard QO is proportional to the extremal 
cross section of the Fermi surface (here directly related to p = 2 ^ 
with the area S = nkp). We find that even if p is inside the gap 
(blue dashed) or in the fiat band region (red) there are well defined 
QO which are directly proportional to the area picked out by the 
intersection of the unhybridized bands (here directly proportional to 

w/hujff 


0.257 where we can replace F ^ Fq = F ^,n = o'j 
to obtain a generalized LK-like form, which has a simple in¬ 
terpretation as a doping and also temperature-dependent ef¬ 
fective mass renormalization. In order to substantiate our un¬ 
expected findings we reproduce all our results in an unbiased 
numerical tight-binding lattice model calculation. 

The model. We consider non-interacting electrons with dis¬ 
persion e(fc) hybridized (strength ^) with a flat band of com¬ 
pletely localized electrons at energy W. The microscopic ori¬ 
gin of such a model is irrelevant for our discussion but the 
Kondo lattice model relevant for heavy fermion systems is ef¬ 
fectively described by such a simple band structure at tem¬ 
peratures well below the Kondo temperature [27-29]. The 
Hamiltonian is simply written as 


H = 


ik) f 
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with the two resulting energy bands E'^{k) = 

^ |e(fc) -\-W ± ^ ~ -F 7^1 separated by a 

hybridization gap 7 and centered around the flat band 
energy W (blue dashed), see Fig. (1) (ii). If p lies within 
the band gap the system is insulating. Once an external 
magnetic field, B = Bz, m switched on (described by a 
vector potential A) the Landau level (LL) structure is easily 
found for a continuum version of our model by replacing 


with a summation over all possible states including all degen¬ 
eracies. We begin with the zero temperature behaviour 

Q(p,T = 0)=iV$ ^ {E±il)-p}. (6) 


We regularize the divergent sum over E-{1) by introducing 
a maximum chemical potential for that lower branch, = 

2 ^ ^max “1“ 1 yj [«max 1] + [w] I which is simply re¬ 
lated to the maximum occupation rimax of the flat band with¬ 
out a field. Here, rimax is defined relative to the filling of a 
dispersive band e(fc) with Fermi energy p = W which de¬ 
fines an occupied area of the BZ S. For our continuum model 
with e{k) = ^ we simply have S = irkp and the relation 
p = 2 ^ straightforwardly generalizes our results to general 
dispersions e(fc) [3]. 

In Fig. I (i) we show the variation of M as a function of 
magnetic field for different chemical potentials (fixed 7 /FF = 
0.05, rimax = 5 and all our findings are independent of the 
cut-off occupation rimax)- For p far above (below) the gap 
there are the usual sharp QO with periodicity directly pro¬ 
portional to the occupied FS volume, see the black dashed (dot 
dashed) curves. For p inside the gap (blue dashed) or inside 
the flat part of the bands (red) we still find well defined anoma¬ 
lous QO of comparable amplitudes. However, now these QO 
have a periodicity hence a BZ area defined by the inter¬ 
section of the unhybridized bands! For larger values of 7 /FF 
(not shown) the amplitude of QO are strongly suppressed for 
smaller magnetic fields but as long as hujc ^ 7 they remain 
observable. 

Effect of temperature. Next, we study the temperature de¬ 
pendence which can by easily calculated for free electrons 
from 0(/i, T = 0) via the convolution [3] 

n{p, T) = - r 0)de (7) 

J-00 0^ 


with the derivative of the Fermi function — = 

2T [ i-r ^ h Strongly peaked at^ = p with a width 

set by temperature. The advantage of this expression is its in¬ 
tuitive interpretation: it is a weighted average over different 
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Figure 2: (color online.) Temperature dependence of the damping 
factor R(T). In (i) it is shown for = 0.7 and for different val¬ 
ues of the chemical potential, ^ = W = Sji, parametrized by ^. 
Dashed lines are calculated from the approximate Ro{T) in which 
r is replaced by Tq = r(n = 0). The inset (ii) shows the same 
for a different value = 0.2. In this case Rq{T) always co¬ 
incides with the exact R{T). Note that for large values of ^ the 
standard Lifshitz-Kosevich behaviour 7?lk (red dashed) is quickly 
approached. 


chemical potentials from a window proportional to tempera¬ 
ture. For standard QO different p correspond to different pe¬ 
riods, hence increasing T always damps the sharp amplitudes 
via dephasing. Evaluating Eq. (6,7) numerically, we find that 
this is not the case for our system: e.g. for jj, = W inside the 
gap we find an initial increase of the amplitudes up to a max¬ 
imum at r Ri 7/4 before damping sets in (not shown). This 
arises because in the temperature average over different p all 
QO have the same periodicity (at least for low T) preventing 
dephasing, however those from regions in the flat part have a 
larger amplitude. 

Eor an analytical calculation of the T-dependence we fol¬ 
low earlier work [12,26] using a finite temperature description 
in terms of Matsubara frequencies a;„ = 2'KiT{n + 7 ). The 
oscillatory part of the grand canonical potential takes the form 
0 (^,T) = nN^J2T=i with l*{n) be¬ 

ing the EL index which defines the pole of the Greens function 
G{iLUn,l) = {iujn — [E±{1) — We write fi = W + 6fi 

and find a single I* to obtain 
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where we have neglected a small n- and 5/r-dependence of the 
real part of I* which only slightly modifies the periodicity but 
not the damping; E is defined below Eq. (2). Now 

differentiating w.r.t. magnetic held and in the limit 
1 we obtain the final result for the first harmonic fc = 1 of the 


Figure 3: (color online.) QO oscillations for a tight-binding lattice 
model with W = —2.3 and 7 = 0.2 (all energies in units of t). In the 
main figure (i) for T = 0 it is confirmed that as long as the chemical 
potential is far away from the gap (black dashed and dot dashed) the 
FS area, 5'*^° extracted from the QO period nicely reproduce the area 
obtained from the relative BZ area of the FS for zero field, see 
legend. If p is inside or close to the gap (green, blue, red dashed) 
we find anomalous QO as before with a periodicity not related to 
In the inset (ii), we extracted the temperature dependence R{T) 
by calculating the difference between a consecutive minimum and 
maximum of M as a function of temperature, which confirms the 
analytical behaviour of Eq. (2), compare to Fig. 2. 


magnetization: 
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with the damping factor R{T) given in Eq. (2). 

In Eig. 2 we plot representative curves of R{T) which fully 
capture the behaviour we have found by numerically evaluat¬ 
ing Eq. (6,7). Eor a chemical potential inside the gap (^ = 0 
black curves) there is an increase of the amplitudes up to a 
maximum T which is set by the energy scale of the hybridiza¬ 
tion 7 itself. The total (relative) height of the maximum in¬ 
creases (decreases) for smaller [see inset (ii)]. Eor larger 
or smaller fillings a characteristic steep increase of the ampli¬ 
tude at a scale T r; 5/r/lO is observed. The simple approx¬ 
imate formula Ro{T), see Eq. (2), in general reproduces the 
behaviour of i?(T) for sufficiently large temperatures [dashed 
curves in (i)]. Eor small values of it fully captures the 
exact result as shown in the inset (ii). 

Lattice model. So far our theory was restricted to a con¬ 
tinuum description, requiring regularization of the flat band 
occupation. To confirm our findings for a microscopic model, 
we have performed a full lattice tight binding calculation. We 
consider a model of spinless electrons on a square lattice, with 
Hamiltonian 

H = - +h.c.) 

{id) 















































4 


The magnetic field is incorporated in the phases of the nearest- 
neighbour hopping parameters tij via the usual Peierls sub¬ 
stitution. These itinerant electrons are coupled locally to a 
second completely localized orbital with on-site energy W at 
each site. The magnetic flux through the magnetic unit cell 
of size LxLy is quantized to multiples of the elementary flux 
quantum <1) = L^LyB = m^Q. We study the system at a 
series of magnetic fields for which Ly = 2 and there is an in¬ 
teger Lx such that the flux $ = <i>o. For each field the Hamil¬ 
tonian is easily diagonalized as before, but now the maximum 
occupation of the flat band is fixed by the total number of 
lattice sites. The QO are directly calculated from the grand 
canonical potential, Eq. (5). 

In Fig. 3 we show the QO of the magnetization which we 
obtain from our lattice simulation. We not only recover the 
anomalous dHvAe at T = 0, see main panel (i), but we also 
confirm the peculiar temperature dependence R{T) of the am¬ 
plitudes, see inset (ii). If the chemical potential lies in the flat 
part of the band such that (5/i <C IE we recover the peculiar 
upturn of the amplitudes towards the lowest T. 

Discussion and conclusion. We have shown that, at odds 
with the canonical understanding of QO in metals, a sim¬ 
ple model of itinerant electrons coupled to a flat band can 
lead to clear QO even in the complete absence of a FS. We 
find strong deviations of the temperature dependence from the 
usual FK theory and derived analytic expressions which can 
be tested in future experiments. We believe that our results 
are most promisingly applicable to certain heavy fermion ma¬ 
terials whose properties well below the Kondo temperature 
are effectively described by a band structure similar to our 
model [27-29]. In that context it is worth pointing out that 
our theory has its most prominent deviations from the FK de¬ 
scription in a regime in which the cyclotron frequency, fiwc, is 
larger than the hybridization strength 7 as well as the activa¬ 
tion gap 7 ^/ 4 IE - a condition fulfilled at least by some heavy 
fermion materials. 

Interestingly, the main features of our peculiar temperature 
dependence were already observed in heavy fermion com¬ 
pounds in two of the rarely available experimental examples 
deviating from FK theory: Amplitudes of some frequencies of 
the dHvAe in CeCoIns display a clear maximum at a nonzero 
temperature of 100 mK, which has been attributed to a fine 
tuned spin-dependent mass enhancement. Most recently, the 
tentative topological Kondo insulator SmBe [30], for which 
the appearance of QO itself despite the opening of an activa¬ 
tion gap [15] (as seen in transport) has been a puzzle, does 
show QO with a very strong increase of intensity below IK 
signaling the presence of a second low energy scale in the 
system [22]. Although, the latter is likely an interaction ef¬ 
fect it is interesting to note that in our non-interacting theory a 
chemical potential not in the gap but just touching one of the 
heavy bands (\Sfi/j\ > 0) [31] sets a new energy scale and 
gives a very similar temperature dependence with a steep in¬ 
crease of the amplitudes at very low temperatures, see Fig. 2. 
For the actual material SmBg it is more likely that our sce¬ 
nario just explains why there are QO in this Kondo insulating 


system at all but the incorporation of self energy effects into 
our theory, which will introduce a second energy scale from 
coupling to collective modes, is a promising route for future 
investigations. In addition, it is an open question for future re¬ 
search, whether certain semiconductors with small direct band 
gaps could also display similar anomalous dHvAes. 

Despite many decades of intense research on the dHvAe we 
have demonstrated that it still holds surprises - there can be 
be QO even in insulating systems. Beyond a mere curiosity 
the interest in standard FK-like QO derives from its capac¬ 
ity of accurately determining FSs. Similarly, we anticipate 
that our anomalous dHvAe applicable to heavy Fermi liquids 
will be useful in the future for determining hybridization gaps 
(proportional to the Kondo coupling) by measuring the tem¬ 
perature of maximum amplitudes. 
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magnetic field on the square lattice. In addition to itinerant 
delocalized electrons with a dispersion e{k) = —2t cos{kx) — 
2t cos{ky) there is a second completely localized orbital with 
on-site energy W at each site. Both d.o.f. are hybridized and 
still described by Eq. (4). The system is diagonalized as be¬ 
fore but now the maximum occupation of the flat band is fixed 
by the total number of lattice sites. A magnetic field corre¬ 
sponding to the discrete vector potential A = Bxy is incorpo¬ 
rated into the tight-binding hopping parameter t via the usual 


i r®’*' A-dr 

Peierl’s substitution e The 

magnetic flux through the magnetic unit cell of size L^Ly 
is quantized to multiples of the elementary flux quantum 
<!> = LxLyB = m^Q. We work in a gauge with fixed Ly = 2 
and m such that we have a one dimensional unit cell of length 
L proportional to B. (We put = L and use an elementary 
flux quantum with units <i)o = 27r.) 

In the following, the upper index I labels the position in¬ 
side each unit cell and (x, y) label the position of the unit 
cells. Itinerant (localized) electrons are created by operators 
y (fx,y)- The lattice Hamiltonian takes the form H = 
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For each field B we have 


a translationally invariant system with unit cells of 
length L. Then, we use a Fourier transform cl. y = 


E ik^x—ikyy d 
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and a corresponding spinor 
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energies are easily found from the resulting 2L x 2L quadratic 
form. 




such that the 


SUPPLEMENTARY MATERIALS 
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Here, we present details of our tight-binding lattice model 
calculation. We use a model of spinless electrons in a uniform with the matrix 
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and the diagonal LxL matrices = ? 1 and Hy/ = W1. potential, Eq. (5). 
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Finally, QO are directly calculated from the grand canonical 











